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Abstract 

It is well established that the mechanical environment influences cell functions in health and disease. Here, we address how 
the mechanical environment influences tumor growth, in particular, the shape of solid tumors. In an in vitro tumor model, 
which isolates mechanical interactions between cancer tumor cells and a hydrogel, we find that tumors grow as ellipsoids, 
resembling the same, oft-reported observation of in vivo tumors. Specifically, an oblate ellipsoidal tumor shape robustly 
occurs when the tumors grow in hydrogels that are stiffer than the tumors, but when they grow in more compliant 
hydrogels they remain closer to spherical in shape. Using large scale, nonlinear elasticity computations we show that the 
oblate ellipsoidal shape minimizes the elastic free energy of the tumor-hydrogel system. Having eliminated a number of 
other candidate explanations, we hypothesize that minimization of the elastic free energy is the reason for predominance of 
the experimentally observed ellipsoidal shape. This result may hold significance for explaining the shape progression of 
early solid tumors in vivo and is an important step in understanding the processes underlying solid tumor growth. 
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Introduction 

Tumorigenesis and solid tumor growth are associated with 
altered mechanics in the tumor's environment such as increased 
matrix stiffness and growth-induced mechanical stress [1-4]. Such 
altered mechanical properties and stresses influence cancer cell 
behaviors such as growth, survival, organization [1—3], and 
tensional homeostasis [4], aid in invasion [5] and effect gene 
expression [6]. Additionally, the tumor's size and morphology may 
be affected [7-9]. 

For example, both in vivo and in vitro, solid tumors are often 
described as being ellipsoidal in shape. In vivo, the growth of solid 
tumors in ellipsoidal shapes is common across tumor classifications 
and tissue of origin [10-16]. This ellipsoidal growth of tumors can 
be emulated in vitro. In a model of the early, prevascular stage of 
tumor growth, tumor cells are embedded and allowed to grow in a 
tissue-mechanics-mimicking hydrogel [7-9]. In such an environ- 
ment, tumor shape has been shown to suffer a loss of symmetry 
from spherical to ellipsoidal for two different tumor cell lines [7,9] 
(murine mammary carcinoma and rat prostate carcinoma, 
respectively). Since the hydrogel is biochemically inert to cellular 
attachment, the effect of the mechanical environment on tumor 
growth is isolated. 

The widespread occurrence of the ellipsoidal tumor shape in 
vivo and the fact that it can be reproduced with multiple cell lines 
in vitro with only mechanical, not biochemical, constraints suggest 
that the mechanical environment strongly influences tumor shape 



development. However, the physical conditions driving ellipsoidal 
tumor growth have not yet been described. Here, we report on our 
study of the mechanics of tumor growth using an in vitro model of 
prevascular tumor growth in chemically inert hydrogels. 

Materials and Methods 

Cell culture and (immuno)fluorescence stainings 

Reagents were from Gibco (Carlsbad) unless otherwise stated. 
Cells — human colon adenocarcinoma (LS174T from ECACC, 
Porton Down) — were propagated in Bio-Whittaker Eagle's min- 
imum essential medium (Lonza, Walkersville) supplemented with 
10% fetal bovine serum, 1% Pen/ Strep, and 1% MEM non- 
essential amino acids. Cells were split prior to becoming confluent 
using 0.25% Trypsin-EDTA and counted with a Z2 Coulter 
Counter (Beckman Coulter). For fluorescence imaging, cell nuclei 
were stained with DAPI (Serva Electrophoresis, Heidelberg) or 
Hoechst (Invitrogen), Alexa Fluor 568 Phalloidin (Invitrogen) was 
used for immunofluorescent staining of the actin, and Alexa Fluor 
594-conjugated E-Cadherin antibody (Cell Signaling Technology). 

Agarose hydrogel preparation, tumor embedding, and 
growth conditions 

Agarose powder (type VII, Sigma, St. Louis) was dissolved in 
heated deionized water. The liquid agarose had twice the agarose 
concentration (in wt./vol.%) than intended for the experiment. 
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Diluting with warm cell culture medium at the time of mixing with 
the cell suspension established the final concentration as well as a 
nutritive environment for the cells. 

Cell-laden gels were maintained in 6-well cell-culture inserts 
with 1-u.m porous membrane bottom (BD-Falcon, Franklin Lakes). 
A cell-free layer of agarose liquid (0.5 mL) was first deposited on 
the bottom of the well insert and allowed to gel before adding the 
cell-containing second layer (3-4 mL). After the second layer of 
agarose had gelled, cell culture medium was added both on top of 
the gel and in the well containing the insert (5-7 mL in total). 

Preproduced tumor spheroids [17,18], formed by the hanging 
drop method [19], were directiy injected with a micropipette tip 
into the gel after it was added to the well insert but before it gelled. 

Tumor Imaging 

After the spheroids were embedded in the gels, their develop- 
ment was monitored for time periods from three weeks up to two 
months by acquiring phase contrast images (Axio Observer.Zl, 
Carl Zeiss AG, Jena), generally at 48-hour intervals. The 
projections of the tumors were imaged in this manner. Only 
tumors that were well away from any boundary of the cell culture 
well were imaged. Since the walls of the well are angled, there is 
optical distortion preventing imaging of tumors growing adjacent 
to the walls. Tumors located within approximately one tumor 
radius away from the top or bottom of the well tended to grow 
with their major axes aligned with the gel boundaries, but the 
experimental setup was not conducive to precise measurements of 
boundary-tumor distances. It is for this reason that only tumors 
well away from the boundaries were measured and considered 
here. 

Confocal fluorescence imaging was performed using a Leica 
TCS SP5 X on the DM 6000 CFS upright microscope (Leica). 

Light sheet fluorescence microscope imaging (for Figure 1) was 
performed on a Zeiss Lightsheet Z. 1 microscope with the help of 
Dr. C. Schwindling at the Zeiss Microscopy Labs in Munich, 
Germany and (for Figure 2) on a custom-built light sheet 
fluorescence microscope in the lab of Prof. H. Schneckenburger 
with the assistance of S. Schickinger at the University of Aalen, 
Germany [20]. 

Finite element computations of tumor growth as a finite 
strain, nonlinear elasticity problem 

The foundation of our computations lies in the definition of the 
deformation gradient tensor, F = l + du/ dX, where u is the 
displacement field vector, and 1 is the isotropic tensor. To model 
the kinematics of growth, we adopt the elasto-growth decompo- 
sition (see Refs [21-23] and [24]) F = F e F g , where the growth 
tensor F g = diag{oi\, 0.2,0-3), and F e is the elastic part of the 
deformation gradient. Here a.\, 0:2 , a 3 ar e growth ratios in the 
Cartesian directions defined by the three axes of a tumor that, in 
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Figure 2. The various projections of oblate ellipsoidal tumors. 

a. Projections (maximum intensity) of two oblate ellipsoidal grown in 
0.5% agarose gel and imaged with a light sheet fluorescence 
microscope from four different directions, b. Rotational sequences of 
the 3D renderings. The 0° orientation is marked with an asterisk in a. 
Complete sequences are available as Movies S3 & S4. Fluorescence 
signal is from Hoechst-stained nuclei. Scale bars are all 50 am. 
doi:1 0.1 371/joumal.pone.01 03245.g002 

general is ellipsoidal in shape. We model the tumor and gel as soft, 
isotropic, nearly incompressible materials, using a neo-Hookean 
strain energy density function: 



<KC e ) = - A(det C e - 1) - - (~ + ^ In det C e + - n(trC e - 3) 

where C e = \ (F eT F e — l) is the right elastic Cauchy-Green 
tensor, and X, ji are elastic constants called Lame parameters, 
which are related to the Young's modulus, E, and Poisson ratio, V, 
by 
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Figure 1. 3D rendering of an oblate ellipsoidal tumor. The 3D 

rendering, created using images taken with a light sheet fluorescence 
microscope, is rotated about the vertical axis in this sequence of 
images. Blue: Hoechst-stained nuclei; Red: E-cadherin. Scale bar is 
90 urn. 

doi:10.1371/journal.pone.0103245.g001 



We have carried out finite element computations of the 
nonlinear, finite strain elasticity problem of the growth of initially 
stress-free spherical and ellipsoidal tumors in a constraining gel 
with the nearly incompressible neo-Hookean strain energy density 
function, using established methods [21-23,25]. The neo-Hoo- 
kean function is derived from statistical mechanical principles 
accounting for the underlying Gaussian network of polymer chains 
[26], which form the gel and the extracellular matrix in the tumor. 
The bulk to shear modulus ratios for tumor and gel were taken to 
be Ktum/t*tum = K gel/Pgel = 5Q' t0 represent the near incompressi- 
bility of both. This corresponds to X tum j fi tum = X ge i j \i gej = 49. Our 
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computations provide the Cauchy stress tensor, ay, and the elastic 
deformation gradient, F|. The latter enables us to compute the 
elastic free energy. 

The initial tumor was modeled as an ellipsoid of semi-axes 
a\, ai, «3, where a\, «2 lay in the range 50-72 um, while 03 was 
determined by requiring all the tumors to have a fixed initial 
volume, imposed by specifying a^aja^ = 1.25 x 10 5 urn . The 
initial ellipsoids thus had aspect ratios of axes varying between 1 
and 3. The ellipsoids were embedded in a cubically shaped gel of 
side 2000 um. The tumor-gel interface allows transmission of 
normal and tangential tractions. This models the effect of bonding 
between the ECM deposited by the tumor cells, and the gel. The 
boundary value problem of nonlinear elasticity was driven by 
specifying the growth tensor to be isotropic, 6(1=6(2=6(3, with 
final growth volume ratio det F g = a\ 6(26(3 = 5 or 10. We specified 
displacement boundary conditions, u n = 0 on five of the six 
boundaries of the gel, while on the sixth we used traction-free 
conditions an = 0. Thus, we modeled the gel in a well with one free 
surface, as in our experiments. We confirmed that all stresses 
decayed by at least two orders of magnitude at the gel boundaries, 
relative to their maxima, thus ensuring that the extent of the gel 
and its shape had no influence on the computations. 

In the finite element computations, the restriction of the tumor- 
gel interface to element edges would have two drawbacks: (a) It 
would require excessive mesh refinement to approximate the 
curved interface, (b) Even with a high degree of mesh refinement, 
stress singularities would arise at the edges and vertices of the 
hexahedral elements that would lie at the tumor-gel interface, 
because of the discontinuities in strain (from tumor growth) and 
elastic constants. Both these limitations can be mitigated if the 
finite element formulation can be extended to allow the tumor-gel 
interface to intersect an element. Such methods are well-known in 
the finite element literature. Our implementation is based on the 
enhanced strain formulation described in Garikipati & Rao (2001, 
[25]), which we have extended to three-dimensional hexahedral 
elements. We note that this method allows the tensors F and F g to 
be discontinuous within elements that contain the tumor-gel 
interface, as dictated by the mathematically exact kinematics of the 
problem. 

We have implemented the finite element formulation oudined 
above in the open source code deal.ii [27]. The stress fields and 
energy surfaces presented as results in this paper were obtained by 
this formulation. Each point on the energy surfaces is the result of 
one computation as oudined above. Typical computations 
involved ~ 300,000 elements and ran for ~30 hours on a node 
with 16 cores, 4GB RAM and a clock speed of 3 GHz. 

Results 

Experimental tumor growth model 

Colon cancer cells (LS174T cell line) were incorporated, either 
as a single-cell suspension or pre-produced tumor spheroids, with 3 
cc of liquid agarose before it gelled. The density of tumors and 
concentration of the agarose gel were controllable. The number of 
tumors in the gel ranged from one to the number that developed 
after inoculation of 10,000 cells/cc (an estimated 50% to 70% of 
embedded "single" cells do not form tumors). Agarose concen- 
trations from 0.3% to 2.0% were used to study tumor growth in 
gels that were in the range of being just more compliant to just 
stiffer (0.3% agarose and 0.5% agarose: 0.3±0.2 kPa and 
0.7 ±0. 1 kPa) or significantly stiffer (1.0% agarose and 2.0% 
agarose: 4.0±0.5 kPa and 24±3 kPa) than pre-produced tumor 
spheroids (0.45 ±0.03 kPa) [28]. All growth experiments were 
conducted in 6-well porous-bottomed cell culture inserts and only 



tumors that were well away from the gel boundaries were 
measured. 

The key parameters measured during growth were tumor size 
and shape. Tumor size, which increases with time after 
embedment [8,9], was affected by both tumor density and agarose 
concentration. Tumor shape, however, was consistently observed 
to develop to oblate ellipsoidal (semi-axes a\ > xa2>a^, Figure 1) 
regardless of tumor density. In the case where pre-produced tumor 
spheroids were embedded that had initial projected tumor aspect 
ratios between about 1.2 and 1.7 (Figure SI), this shape formed 
within one week after embedment (Figure SI; Movies SI, S2). We 
confirmed that it did not form due to collapse of a necrosed core 
(Figure S2). The only parameter affecting this shape development 
was the agarose elastic modulus: when the agarose gel was stiffer 
than the pre-produced tumors, the latter grew as oblate ellipsoids. 
When the gel's elastic modulus was below that of the pre-produced 
tumors, approximately spherical tumors grew, with diffuse 
boundaries (Figure S3). 

Although characterization is presented here for LS174T colon 
adenocarcinoma cells, we also have observed the oblate ellipsoidal 
shape of tumor development with HeLa cervix adenocarcinoma 
cells embedded in 0.5% agarose hydrogels (Figure S4). 

Tumor shape and orientation characterization 

The three-dimensional shape of the tumors was investigated 
using light sheet fluorescence microscopy (Figures 1 and 2) and 
physical sectioning of gels to obtain at least two perpendicular 
views of individual tumors with a confocal microscope (Figure S5). 
Full 3D measurements were made of 1 1 tumors in 0.5% or 1.0% 
agarose (Table SI in File SI). Eight different null hypotheses were 
devised to test whether these shapes were oblate ellipsoidal and the 
Bonferroni adjustment was used to correct for the increased risk of 
type I errors (false positives, further details in File SI). The sum of 
the p values for all 8 tests was 0.0064 (Table S2 in File SI), well 
under the set statistical significance level for each test (here, p/ 
8 = 0.00625, obtained for p = 0.05). Furthermore, the average 3D 
aspect ratio («i/«3) of tumors in the 1% agarose gel was 2.7±0.3, 
confirming uniformity in oblate tumor shape. 

Three-dimensional reconstructions (Figure 2b; Movies S3, S4) 
of light sheet-imaged tumors revealed that the oblate ellipsoidal 
tumor shape was correlated with a wide range of projected 
elliptical shapes, from a circle [ci\ x « 2 ) to an ellipse of the 
maximum achievable aspect ratio (ai/a^), with a range of smaller 
aspect ratios in between (Figure 2a). The length of the major axis 
of the ellipse that is projected by an oblate ellipsoid will always be 
equal to the length of the largest axis of the oblate ellipsoid itself 
(Figure S6 in File SI, the outline of a proof is presented in File SI). 
Therefore, all possible projected aspect ratios of such an oblate 
ellipsoid may be calculated by rotating it around one of its major 
axes (Figure 3a). 

Since the oblate tumor shape is robust and its maximum 3D 
aspect ratio is relatively uniform for tumors grown in 1 % agarose 
gel, we have used the projected aspect ratio as a measure of tumor 
orientation — defined by rotation around one of its major axes 
(Figure 3a) — for many tumors. In several experiments with a cell 
inoculation density of 2500 cells/cc in 1% agarose gel, a 
maximum aspect ratio of 3 was observed after one month of 
growth, which is consistent with the average 3D aspect ratio 
measurement of 2.7±0.3 under the same conditions. The 
distribution of projected aspect ratios, and therefore tumor 
rotations, is shown in Figure 3b for one experiment. After 
removing the gel and sectioning it to image the tumors from the 
side, a similar distribution of projected aspect ratios resulted 
(Figure 3b). This lack of preferential tumor orientation indicates 
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Figure 3. Oblate ellipsoidal tumors have a similar distribution of orientations regardless of viewing direction, a. The relationship 
between tumor rotation angle, a, and projected aspect ratio, a\ja p ,j, as an oblate ellipsoid with maximum aspect ratio of 3 (ai/ai = 3) is rotated 
about its a\ (,yi) axis. b. The tumor rotation angle was calculated from the projected aspect ratio (part a) for tumors grown in 1% agarose for 30 days. 
Top image: observation perpendicular to the plane of the cell-culture well, Bottom image: side view perpendicular to a physical cross-section of the 
gel made with a scalpel blade (out-of-focus cut marks are visible in the gel) Scale bars are 200 |im. 
doi:1 0.1 371 /journal. pone.0103245.g003 



there is no influence of an externally applied field, such as a 
mechanical stress field, which could bias all the tumor ellipsoidal 
axes to the same direction. 

Nonlinear elasticity finite element computations 

If an elastic field is associated with the ellipsoidal shapes, such a 
field must bear a random relation to position since the major axes 
of the ellipsoidal tumors are randomly oriented with respect to the 
tumor's position in the gel (Figure 3b). One such example is a 
tumor growth-induced stress field, which arises in the tumor and 
gel. If it were unconstrained by the gel, uniform growth would 
cause a stress-free expansion of the tumor, with possible shape 
changes, such as of a sphere to an ellipsoid. It is represented 
exactly by the finite growth tensor F g = diag{a\, 0.2,0.3). The 
determinant gives the growth volume ratio: det F g = ai«2«3. 
When constrained by the gel, however, the actual expansion of the 
tumor differs from F g , and is given by the deformation gradient 
tensor, F. It satisfies the previously introduced elasto-growth 
decomposition F = F L 'F g [21] where detF<detF 8 , and F e 
causes the growth stress. The random position, orientation, and 
growth tensors of the tumors make the growth stress tensor due to 
each tumor also random in magnitude and orientation. 

The stresses in the tumor and gel are reported in Figure 4 for a 
spherical tumor and an oblate ellipsoidal tumor 
[a\ =cii, a\/a^ = 3), both growing isotropically (ai=«2=«3) in 

gels that are ten times stiffer (fltum / Pgel = V 1 "), matching our 
experiments for tumors in 1 % agarose gel. Growth leads to a 
compressive stress inside the tumor when it is constrained by the 
gel. For a spherical tumor growing isotropically, the compressive 
(negative) stress is uniform in every direction 
(a 1 1 = ff22 = 033 = Osph)' However, for an ellipsoidal tumor growing 
isotropically, the stress components within the tumor are no longer 
equal: For an oblate ellipsoidal tumor [a\ =fl2> a l/ a 3 > 1), the 
compressive stresses along the a\ and a-i axes satisfy (J\ \ = G22, and 
are greater than that along the short axis, |cri 1 1 > | (T33 1 . For a 
growth tensor F g , these compressive stress magnitudes can be 
computed to be continuous at the tumor-gel boundary, and to 
decay outside the tumor. Tensile normal stress components are 
induced in the hydrogel tangential to the tumor-gel boundary, but 
no cracking of the hydrogel was detected as a result (Figure S7 in 
File SI, for further details see File SI). 

We computed the elastic free energy, of the tumor-gel system 
for large isotropic growth of the form 7 rg = <i/ag(5 1 < /3 ,5 1 < /3 ,5 1 < /3 ), 



detF ? = 5. A total of 576 boundary value problems were solved, 
each with different ratios a\/a^ and 02/03. The initial ellipsoidal 
volume was fixed by imposing a\Ci2Ci3 = 1.25 x 10 5 um\ Results 
have been plotted in Figure 5 for the elastic free energy versus 
fll/fl3 and «2/fl3 in the interval 1 <fli/fl3, fli/fl3 < 3. Matching 
our experiments with tumors grown in 1 % agarose hydrogel, 
Figure 5a summarizes the elastic free energy results for 

fi tum jli gel = 1/10. The free energy is highest in the neighborhood 

of approximately spherical tumors, 01/123 ~ 1—1.5 and fl2/ fl 3 ~1— 
1.5, and lowest in the limit of oblate ellipsoidal tumors: 01/03, 
02/03 ~3 — the shape most often observed in our experiments. 
The energy surface in Figure 5a also shows that for a fixed value of 
ai/03 (02/03) the energy is minimized for the largest 02/03 (fll/ a 3), 
although this variation is more gradual. This result suggests that, if 

fmm I Pgel < l s the reason for the predominance of the oblate 




-1400 I ' ' ' ■ 1 

0 1 2 3 4 5 

Distance from ellipsoid origin (x 10 2 ), \im 

Figure 4. The stress field created by tumor growth when 

Mtum/Pgel = 1/10. The growth volume ratio det F g = 10. All normal stress 

components, tr„ ( ), are equal in a spherical tumor of initial radius 

50 |j,m. The maximum compressive stress is —1300 Pa. However, in an 

ellipsoidal tumor with axes a\ =cin = 3^3 = 72 |rm, the dn ( ) 

component, has a maximum compressive stress of —1200 Pa, and the 
(733 ( ) component has a maximum compressive stress of - 
1 1 80 Pa. 

doi:1 0.1 371/journal.pone.01 03245.g004 
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Figure 5. Elastic free energy landscapes of the tumor-gel 
system. The landscapes are plotted versus the tumor's ellipsoidal axes 
ratios 02/03 and 01/03. The origin has been shifted, so that the 
maximum energy is zero in each case: (a) For fi lum j fi gL ,i = \/\0. The 

oblate shapes fa/a? = 01/03 = 3) are the low energy states. This surface 
has been generated from 576 separate nonlinear elasticity computa- 
tions of tumors of aspect ratios 01/03 and 02/03 varying between 1 and 

3, growing in a gel. (b) The same as (a), but for Ptum/ftgei = N°te the 

flatness of the landscape relative to (a). Spherical shapes are not 

penalized strongly for Mtum/ Mgel — ^' DUt are strongly penalized for 

/'/™,A'«w =1 / 10 - 

doi:1 0.1 371 /journal.pone.01 03245.g005 

ellipsoidal tumor shape could be associated with the fact that this 
shape minimizes the elastic free energy. 

In vivo soft tissue tumors, however, are observed to be stiffer 
than the surrounding tissue [2,29,30]. While we have not been 
able to replicate this regime in our experiments, we approach it in 
the case of 0.3% agarose gels, which corresponds to 

Ptumj f l gel < ~1- A- s shown in Figure 5b, the elastic free energy 

landscape is much flatter than in Figure 5a, because the energy 
differences are much lower. This suggests that configurations 
closer to spherical are subject to a less stringent free energy 
penalty. (Note: These results — both for the stress field and the 
elastic free energy — remain valid even if linearized, infinitesimal 
strain elasticity is used following Mura's [31] treatment of 
Eshelby's inclusion problem [32,33]. While the formulas are more 
tractable and the results may be directiy plotted without the need 
for finite element computations, the validity of the calculations is 
doubtful because of the finite strains implied by large tumor 
growth ratios.) This finding corroborates the observation in our 
experiments in this case for which the tumors remain closer to 
spherical (Figure S3). 

Discussion 

Monitoring displacements in the hydrogel surrounding ellipsoi- 
dal tumors using co-embedded fluorescent micro-beads, Cheng 
and co-workers [7] sought to explain the observed symmetry- 
breaking that leads to ellipsoidal tumors grown from murine 
mammary carcinoma cells. The strain fields computed from the 
micro-bead displacements were interpreted as showing the 
compressive stress to be greater along the minor axis of the 



ellipsoid. By correlating these fields with tumor shape and 
Caspase-3 activity, the authors concluded that mechanical stress 
was causing a higher fraction of cell death along the minor axis 
and driving the tumor to grow in the corresponding ellipsoidal 
shape. The authors did not, however, identify the origin of the 
stress in their system. 

We were unable to determine any correlation between tumor 
orientation and position in the gel (Figure 3b). As we have 
observed above, if the ellipsoidal shape were associated with an 
elastic field, that field could not be correlated with position in the 
gel because of the observed random distribution of ellipsoidal 
tumor axis orientations. We also have noted that the stress field 
caused by the growth (expansion) of the tumor in the gel is an 
example of such an elastic field. Nonlinear elasticity, however, runs 
contrary to the conclusions of Cheng and co-workers [7]. To 
demonstrate this, we have carried out nonlinear elasticity 
computations to show that the oblate ellipsoidal growth of a 
tumor induces compressive stress components along its major axes 
that are greater in magnitude than the component along the minor 
axis when measured just outside the tumor (Figure 4). This result is 
also obtained with linearized, infinitesimal elasticity. If compres- 
sive stress suppresses growth by signaling different cell growth 
and/or death rates, the major axes ought to suffer this suppression 
more sharply than the minor axis. Such a sequence of events 
would favor the spherical shape, which clearly does not happen in 
our experimental studies. The experimental evidence suggests that 
the growth-induced compressive stresses, which are higher along 
the major axes, do not suppress the growth of ellipsoidal tumors. 
Therefore, such shapes must be controlled by another quantity. 

Over the course of our studies, we computationally tested 
several mechanisms that could drive tumor growth into an oblate 
ellipsoidal shape. These included (a) growth along a compliant 
layer, (b) the rise of a highly proliferative subpopulation of cells, (c) 
stress-driven migration of cells, and (d) elastic free energy 
minimization. Through extensive tumor growth modeling based 
on our past work [21-23] and that of Casciari et al. [34], we 
concluded that the mechanisms in (a)-(c) could account for small 
perturbations of the tumor shape toward an ellipsoid, but not for 
the strongly ellipsoidal shapes that we have observed (see reference 
[1-4,23] in this regard). The study of the elastic free energy of the 
tumor-gel system, however, has proved more promising. Indeed, 
for the case of a more compliant tumor in a stiff matrix 

{P-tumj P l gel = ^l^) a steep energy surface results in which the 
elastic free energy is significantly minimized for oblate ellipsoidal 
tumor shapes (Figure 5a). Although the oblate ellipsoidal shape is 
still the energy-minimizing shape if the tumor and gel have equal 
shear moduli, the free energy penalty for remaining closer to a 
sphere is just about a fifth of the penalty for the case 

/*tum/f t gel= l / W ' 

We confirmed that the orientation of the semi-axes of the 
ellipsoidally growing tumor had no influence on the total strain 
energy of the tumor-gel system, as long as the gel boundaries were 
sufficiently distant. We also confirmed that the location of the 
tumor relative to the gel boundaries had negligible effect on the 
elastic free energy. For a tumor of a given shape the elastic free 
energy varied by less than 1% when its position in the gel was 
varied. This included positions that placed the tumor edge within 
one radius, or one semi-major axis of the gel boundary. This result 
is expected also from linearized elasticity. Given these computa- 
tional results, we suggest that any tendency for tumors to grow 
parallel to gel boundaries is due to mechanical contact resulting in 
a force that deforms the growing tumor into such a shape. We did 
not report contact mechanics computations of the elastic free 
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energy for tumors growing in contact with the gel boundary. The 
elastic free energy in such contact mechanics computations is very 
sensitive to the gap function chosen to impose contact. As noted 
above, our experimental set up was not conducive to precise 
measurement of tumor proximity to gel boundaries, because of 
optical distortion at the boundaries. We therefore are not in a 
position to make a rigorous report on the boundary effect on the 
basis of our experiments. However, as also stated above, all tumors 
reported in our experiments were far from the gel boundaries. 

We note that the qualitative trends observed with the 
formulation of tumor growth as a problem of finite strain 
nonlinear elasticity are reproduced with infinitesimal, linearized 
elasticity following Refs [1-3,31-33]. However such calculations 
using the linearized theory lack the physical consistency and 
quantitative accuracy of the full nonlinear theory. 

Extrapolating the above results to in vivo observations, we 
propose that the observed tendency toward ellipsoidal solid 
tumors, as seen in multiple tissue types [10-13,35], results from 
the lower elastic free energy of the ellipsoidal shape in the early 
prevascular stage of growth when they are compliant relative to 
the surrounding tissue. (Note: While we are not aware of any 
observations of pre-vascular tumors in vivo, we note that the pre- 
vascular tumors in our study have moduli in the range of 0.3- 
0.7 kPa, which is at the lower extreme of soft tissues as reported by 
Yu et al. [2,5].) Subsequently, as these tumors grow larger and 
vascularize, and their microenvironment becomes more fibrotic, 
the elastic modulus of the tumor tissue may exceed that of the 
surrounding tissue [2,6]. In this regime of tumor versus 
surrounding tissue stiffness, spherical tumors are not subject to 
as stringent an energy penalty, and therefore are more likely to be 
found. However, with early prevascular tumors being more 
compliant than their microenvironment, the penalization of the 
spherical shape at this stage ensures that it is the ellipsoidal tumors, 
especially the oblate ellipsoidal ones, that are favored. Our 
nonlinear elasticity computations have further shown that there is 
a stringent energy penalty against shape changes that could be 
effected by non-diagonal growth tensors: F g = diag(u.\,a2,ct-i), 
where at least one of oci, Cf-2 and a 3 is different from the others. 
Because of these energy-dependent mechanisms, ellipsoidal tumors 
are favored over other configurations, and are observed more 
often, across all tumor-gel stiffness ratios. 

Our findings do not suggest what mechanisms of growth could 
actually suppress the high-energy shapes. Indeed, our methods, 
experimental and theoretical, are not designed to answer this 
question. One possibility is that, for tumor-gel systems with a 
greater elastic free energy, the mechanical work performed by a 
growing tumor to store this energy drains the biochemical free 
energy of the cells in the early stages of tumor development. This 
hypothesis would have to be tested against the complexities of 
cancer metabolism [7-9,36] before reaching a firmer conclusion. 
Because agarose gels are biochemically inert, mechanical interac- 
tions play a more dominant role in determining the shapes of 
tumors in our study, by design. In vivo, of course, there are other 
biochemical drivers of tumor growth. However, the results of this 
study would suggest that minimization of elastic free energy is an 
important driver of tumor shape and, due to significant reporting 
of ellipsoidal tumors, possibly plays a significant role in vivo also. 

The importance of this work may prove to be that tumor shape 
plays a role in determining the potential for malignancy. Once it 
has a well-developed necrotic core, a tumor can be modeled as an 
internally pressurized thin shell carrying meridional and azimuthal 
stress resultants in its wall, which is composed of viable cells. These 
stress resultants are known from growth models to control shape 
instabilities of the tumor wall [10-16,37-39], leading ultimately to 
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its breakdown. Thin shell theory [7-9,40] shows that, while 
spherical shells have equal azimuthal and meridional stresses, both 
being uniform over the shell, ellipsoidal shells have non-uniform 
fields with extrema at the poles of the ellipsoid. These points on an 
ellipsoidal tumor could therefore be critical sites for the surface 
instability, breakdown of the wall, and potential cell escape leading 
to malignancy. 

Supporting Information 

Figure SI Individual tumor development. To determine 
the time frame over which the oblate ellipsoidal shape developed, 
individually embedded, pre-produced tumor spheroids were 
followed as they grew in the agarose gel. In such an experiment, 
it was determined that within one week after embedment the 
ellipsoidal shape was detectable and the orientation of the tumor in 
the gel was fixed. Shown here is the time course development of 
the projected tumor aspect ratio of 1 1 different individual tumors 
(sparsely embedded). Diamond symbols mark the aspect-ratio 
measurements for tumors whose final orientation projected narrow 
ellipses, solid circles for tumors whose final orientation projected 
near-circles, and open circles for tumors whose final orientation 
projected wider ellipses. Phase-contrast images show the projec- 
tions from which the aspect ratios were measured for the specific 
tumors and time points circled in grey on the plot. Scale bars are 
500 urn. 
(TIFF) 

Figure S2 In the ellipsoidal cross-section of a three-day- 
old tumor, a necrotic core is not detectable. It was 

confirmed that the tumors took on the oblate ellipsoidal shape 
before the observation of a central necrotic core, ruling out the 
possibility that this shape formed due to mechanical collapse of a 
necrosed core. Three-dimensional projection of 13 image slices, 
taken 2 urn apart in the mid-section of an oblate ellipsoidal tumor. 
The section is shown here tilted about the vertical axis of the 
image by 10 degrees. Nuclei (blue) and actin network (red)of 
LS174T cells embedded in 1% agarose hydrogel. Scale bar is 
30 um. 
(TIF) 

Figure S3 Tumors grown in 0.3% agarose. Left Column: 

Phase contrast images of 18-day-old tumors taken with a 5x 
objective. Scale bar is 500 um. Right Column: Projections at 
two perpendicular orientations of a 15-day-old tumor. The 
fluorescent signal is from the Hoechst-stained nuclei. The scale 
bar is 50 um. The tumor shape in the softest gels is oftenmore 
diffuse and approximately spherical. 
(TIFF) 

Figure S4 Phase-contrast images of the projections of 
four different tumors grown for 9 days from HeLa cervix 
adenocarcinoma cells embedded in 0.5% agarose hy- 
drogels. Elliptical and circular cross-sections, similar to those 
seen from tumors grown from LS174T cells, are observed. The 
scale bar is 100 um. 
(TIF) 

Figure S5 Perpendicular views of 5 tumors grown in 1 % 
agarose. Due to scattering regular confocal microscopy is not 
effective at capturing entire 3D tumor shapes. However, by 
carefully sectioning the gel in the region of a tumor, two 
approximately perpendicular views may be imaged of the tumors. 
Here are some examples. The fluorescent signal is provided by 
Hoechst-stained nuclei. The scale bars are 100 |jm long. 
(TIFF) 
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File SI The Supporting Information file contains the 
supporting text for the statistical testing on measured 
3D tumor dimensions, the oudine of the proof for the 
projection of an oblate ellipsoid always containing its 
major axes, and further details on our observation that 
the hydrogel did not crack due to tensile normal stress 
components induced tangential to the tumor-gel bound- 
ary. The Supporting Information file also contains Figures S6 and 
S7 and Tables SI and S2. 
(DOCX) 

Movie SI Time-lapse progression of an individually 
embedded tumor that developed an elliptical cross- 
section. Scale bar is 500 um. 

(AVI) 

Movie S2 Time-lapse progression of and individually 
embedded tumor that developed an approximately 
circular cross-section. Scale bar is 500 urn. 

(AVI) 
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